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SUMMARY 


The  intended  purpose  of  this  investigation  was  to  perform  a thorough 
evaluation  of  the  CTASS  (Controlled  Towed  Array  Sonar  System)  concept. 

This  evaluation  was  to  be  made  under  expected  typical  operation  con- 
ditions with  respect  to  several  performance  criteria:  detection,  reso- 

lution, and  resistance  to  acoustic  interference.  Because  of  mathe- 
matical complexities  in  the  analysis  of  multiplicative  processing  sys- 
tems, the  evaluation  was  to  be  made  largely  by  simulation.  The  project 
consisted  of  several  stages: 

1.  Collection  of  sea  noise  data. 

2.  Mathematical  modeling  of  CTASS  operating  on  simulated  signal  and 
measured  noise  data. 

3.  Filtering  and  digitizing  of  recorded  sea  noise. 

4.  Numerical  CTASS  performance  evaluation. 

The  first  three  stages  were  fully  completed.  In  the  fourth  stage,  num- 
erical CTASS  evaluation,  the  first  performance  criterion  to  be  evaluated 

was  detection,  as  measured  by  the  performance  index  of  the  output  signal- 
e 

to-noise  ratio  (jj)0*  (For  long  integration  times,  detection  probability 

is  a function  of  (S/N)Q.)  Earlier  analysis  (which  had  been  possible  only 

to  a limited  extent)  had  indicated  that  there  was  come  reason  to  believe 
c 

that  (jj)Q  might  not  show  any  improvement  (indeed,  might  show  a decline) 
as  the  size  of  the  array  increased.  It  was  decided  to  cease  further  CTASS 
evaluation  with  respect  to  the  other  criteria  if  signal  detectability  was 
found  not  to  increase  with  array  size. 

Numerical  CTASS  evaluation  showed  a definite  (and  discouraging)  decrease 

s 

in  (^)  as  the  array  size  increased.  This  was  found  to  be  true  even  at  high 
input°signal-to-noise  ratios.  This  prompted  a further  analysis  using  noise- 
free  conditions  with  an  idealized  sinusoidal  input  signal,  but  assuming 
realistic  integration  times.  This  analysis  showed  a definite  decrease  in 

C 

(-)  as  the  array  size  is  increased,  resulting  in  the  conclusion  that  CTASS 
in°its  presently  conceived  configuration  is  not  a promising  array  process- 
ing system.  All  further  evaluation  work  was  therefore  halted.  Possibly 
some  other  configuration  using  gradient  and  omnidirectional  hydrophones 
would  result  in  better  detection  performance.  A modification  of  the  pres- 
ent simulation  model  could  be  used  in  that  evaluation. 
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INTRODUCTION 


Controlled  Towed  Array  Surveillance  System  (CTASS)  refers  to  a proposed 
passive,  Towed,  Acoustic  Surveillance  System  that  would  employ  a multi- 
plicative array  of  gradient  hydrophones,  as  compared  with  the  conventional 
linear  summed  array  of  omnihydrophones.  CTASS  would  employ  two  arrays  of 
nine  electrodynamic  gradient  hydrophones  and  one  omnihydrophone.  The  two 
arrays  would  form  an  orthogonal  set  housed  in  a small  controllable  hydro- 
dynamic  body  about  10  inches  in  diameter  and  9 feet  long  (figure  1).  Ar- 
ray gain  would  be  obtained  by  uniquely  processing  the  output  of  these  hy- 
drophones . 
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Figure  1.  CTASS  Towed  Vehicle 
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The  signal  processing  would  be  multiplicative  and  enhanced  by  the  inter- 
relationship of  the  placement  and  directivity  of  each  hydrophone.  Be- 
cause of  the  acoustic  characteristics  of  each  hydrophone,  and  the  type 
of  signal  processing  used,  the  system  theoretically  would  achieve  a 
three-dimensional  (pencil-shaped)  beam  pattern  that  would  be  invariant 
over  the  frequency  range  of  10  Hz  to  10  kHz. 

Each  CTASS  array  would  consist  of  an  omnidirectional  hydrophone  and 
2q  + 1 gradient  (or  dipole)  hydrophones.  The  omnielement  and  one  dipole 
would  be  combined  to  form  a cardioid  beam  pattern.  The  remaining  2q  di- 
poles would  be  canted  relative  to  the  cardioid  by  progressively  larger 
angles  (n)>  the  first  pair  or  section  at  + n°.  the  second  section  at  + 

2n°,  etc.  The  array  output  would  be  obtained  by  forming  the  product  of 
the  cardioid  output  and  the  outputs  of  all  2q  + 1 dipoles  (i.e.,  the  out- 
put would  be  a product  of  2q  + 2 factors) . The  element  outputs  would  be 
filtered  with  narrow  band  filters,  and  the  product  would  be  low  pass  fil- 
tered, figure  2.  Features  such  as  simultaneous  processing  of  multiple 
frequency  channels,  simultaneous  beamforming,  and  pencil  beams  would  be 
available. 

The  formation  of  the  beam  may  be  described  briefly  as  follows:  The  first 

dipole  would  have  nulls  at  + 90°  from  the  beam  axis.  The  next  pair  would 
have  nulls  at  + 90°  + n°  and  + 90c  - n°.  The  product  of  these  three  would 
have  nulls  at  90°  - n°,  90°,  90°  + n°»  -90°  - n°,  -90°,  and.  -90°  + n°. 

The  next  pair  of  dipoles  would  have  nulls  at  + 90C  + 2n°,  and  inclusion  of 
these  in  the  product  would  produce  a pattern  with  5 nulls  on  each  side  at 
spacing  n°  apart.  Continuing  in  this  manner,  one  would  obtain  a pattern 
with  2q  + 1 nulls  on  each  side  with  quite  low  sidelobes  between  the  nulls. 
The  pattern  thus  far  would  be  symmetric  fore  and  aft.  The  cardioid  pat- 
tern would  have  a single  null  at  180°.  Multiplication  of  the  product  of 
the  2q  + 1 dipole  outputs  by  the  cardioid  output  theoretically  would  de 
stroy  the  back  lobe,  and  produce  a pattern  with  a single  main  lobe  and 
4q  + 3 nulls. 
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Figure  2.  "N"  Section  Array  Signal  Processing  Block  Diagram 


The  CTASS  method  of  beamforming  is  illustrated  in  figure  3.  For  example, 
consider  a two-section  array  and  a canted  angle.  The  two-section  array 


Figure  3.  CTASS  Beam  Forming  Example 

would  use  five  hydrophones  (four  gradient  and  one  cmr.i)  to  beamform.  To 
form  a beam  in  the  0°  direction,  the  basic  processing  unit  (BPU)  would  be 
formed  with  A and  A'  hydrophones  and  augmented  with  the  outputs  of  hydro- 
phones B,  B',  C,  and  C'.  A beam  could  be  formed  in  the  180°  direction  by 
using  the  same  hydrophones  and  reversing  the  polarity  of  the  output  from 

the  A hydrophone.  Reversing  the  polarity  of  the  BPU's  gradient  hydro- 

c o 

hone  would  shift  the  major  response  axis  from  0 to  18C  . 

A beam  may  be  formed  in  the  n°  direction  by  using  the  hydrophone  whose 
main  response  axis  is  in  the  n°  direction,  to  form  the  BPU;  i.e.,  the  BPU 
would  consist  of  hydrophones  B and  A' . A two-section  array  would  now  be 
obtained  by  augmenting  the  BPU  with  hydrophones  C,  A,  D,  and  B'.  A beam 
could  be  formed  in  the  (n  = 180°)  direction  by  using  the  same  hydrophones, 
but  reversing  the  polarity  of  hydrophone  B. 
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The  net  effect  of  CTASS  beamforming  is  a pattern  that  would  be  invariant 
with  respect  to  the  steered  direction,  and  therefore  should  result  in  an 
array  gain  that  would  be  constant  with  respect  to  the  steered  direction. 

The  single  target  beam  response  of  this  array  showed  that  apparently 
quite  narrow  response  patterns  could  be  achieved  while  having  no  require- 
ment on  element  spacing  (except  in  angle) . CTASS  thus  appeared  to  offer 
something  akin  to  super  gain.  Figure  4 shows  the  comparative  single  target 
characteristics  for  arrays  consisting  of  up  to  six  hydrophone  sections. 


The  signal  processing  properties  of  the  CTASS  system  were  analyzed  under 
several  operational  conditions  (reference  1).  As  far  as  possible,  general 
expressions  were  developed  for  the  output  mean,  variance,  and  signal-to- 
noise  ratio  at  the  output  of  the  averager.  In  addition,  a method  was  de- 
scribed by  which  false-alarm  probabilities  could  be  evaluated  under  cer- 
tain conditions.  Because  of  the  exceedingly  complex  nature  of  the  mathe- 
matical expressions  involved  (theoretical  analysis  of  multiplicative  pro- 
cessing is  very  difficult,  and  the  difficulty  increases  rapidly  with  the 
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number  of  multiplications) , numerical  values  for  the  output  variance  and 
output  signal-to-noise  ratio  were  obtained  only  for  a CTASS  array  con- 
taining two  and  four  elements.  The  results  of  this  analysis  showed  that 
the  output  signal-to-noise  ratio  of  the  CTASS  processor  would  be  superior 
to  a conventional  linear  processor  using  the  same  number  of  sensors.  In 
addition,  the  CTASS  resolution  was  investigated  by  analyzing  array  output 
for  various  operational  conditions  when  two  targets  are  present.  This 
analysis  indicated  poor  resolving  capabilities  when  targets  ware  identi- 
cal (both  in  power  level  and  frequency  content)  and  separated  less  than 
90°.  However,  this  analysis  was  made  on  a broadband  basis,  and  it  was 
felt  that  a narrowband  processing  system  (as  proposed  for  the  operational 
system)  would  provide  the  expected  good  performance. 

The  Applied  Physics  Laboratory  (APL)  of  Johns  Hopkins  University  (refer- 
ence 2) , has  made  a limited  evaluation  of  CTASS,  using  simulated  isotro- 
pic noise  and  discrete  samples  of  target  signals.  APL's  simulation  re- 
sults showed  the  same  discouraging  trend  of  decreasing  detection  prob- 
ability with  increasing  CTASS  array  size.  However,  Electric  Boat  had 
serious  doubts  about  the  validity  of  the  APL  simulation  model  and  the 
method  of  performance  evaluation,  so  the  task  reported  here  was  proposed. 

To  resolve  these  issues,  it  was  proposed  that  a suitable  test  be  evolved 
to  establish  the  concept's  feasibility  in  the  static  mode.  Since  the 
difficulty  of  analysis  arises  from  the  mathematical  complexities  of  the 
noise  terms,  it  was  suggested  that  an  experimental  investigation  be  under- 
taken in  which  real  background  noise  data  would  be  employed.  Since  it  is 
anticipated  that  the  CTASS  would  be  ambient  noise  limited,  the  required 
background  noise  would  consist  only  of  sea  state  noise.  The  recorded 
noise  data  would  then  be  used  in  a computer  simulation  that  would  para- 
metrically evaluate  the  CTASS 's  usability  for  antisubmarine  warfare. 

This  report  describes  that  investigation  and  includes  descriptions  of  the 
test  program,  the  simulation  model,  and  the  results. 

OBJECTIVE 

« 

The  objective  of  this  project  was  to  demonstrate  the  feasibility  cf  the 
CTASS  concept  as  an  alternative  to  conventional  towed  sonar  arrays. 

To  that  end,  a two-stage  static  mode  test  was  undertaken.  In  the  first 
stage,  a horizontal  CTASS  array  was  used  to  acquire  real  background  noise 
at  the  Exhuma  Sound  Acoustic  Range.  During  the  second  stage,  data  re- 
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cordings  of  individual  hydrophone  receptions  were  used  to  evaluate  the 
CTASS  array's  ability  to  detect  targets  through  computer  simulations 
carried  out  at  Electric  Boat  Division's  Data  Reduction  Center.  Beam 
patterns,  beam  forming  capabilities,  and  array  sensitivities,  as  affected 
by  array  sections  fi.e.,  number  of  hydrophones),  were  to  be  evaluated 
and  compared  with  those  obtained  from  a linear  array  of  an  equal  number 
of  omnielements  when  narrowband  target (s)  signals  are  present. 

TEST  PROGRAM 
Test  Array 

The  array  configuration  as  tested  appears  in  figure  5.  This  array  consis- 
ted of  nine  gradient  and  six  omnidirectional  hydrophones.  These  hydrophones 


Figure  5.  CTASS  Static  Test  Array 


were  components  of  hydrophone  assemblies  from  dismantled  DIFAR  buoys. 

The  orientation  of  the  hydrophones  within  the  array  is  highlighted  in 
figure  6,  which  indicates  that  the  angle  used  is  10°,  and  that  the 
hydrophone  elements  are  mounted  at  6-in.  intervals.  The  array  has  an  ap- 
propriate number  of  hydrophones  to  evaluate  its  performance,  with  as  many 
as  four  hydrophone  sections  for  a broadside  beam  (nine  gradient  and  one 
omnihydrophones) , three  hydrophone  sections  for  a beam  at  + 10°  (seven 
gradient  and  one  omni hydrophones) , and  two  sections  for  a beam  at  + 20° 
(five  gradient  and  one  omriihydrophor.e)  . All  hydrophones  were  calibrated 
for  sensitivity  in  the  100  to  2400  Hz  band.  In  addition,  beam  patterns 
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were  measured  for  all  of  the  gradient  hydrophones  at  150,  300,  and  2000  Hz. 
Figure  7 is  a typical  calibration  beam  pattern  for  a gradient  hydrophone. 
Typical  sensitivities,  at  the  indicated  frequencies  for  each  of  the 
hydrophones,  appear  in  figure  6.  Each  hydrophone  assembly,  prior  to 
calibration,  was  pressure  tested  to  a pressure  equivalent  to  a 200-ft 
submergence . 
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Figure  6.  CTASS  Test  Array  Hydrophone  Orientation  and  Sensitivity 


Data  Acquisition 

The  CTASS  static  test  array  was  used  to  measure  ambient  sea  noise  at  the 
Exhuma  Sound  Acoustic  Range  during  the  week  beginning  28  March  1976.  The 
array  was  suspended  100  ft  below  the  sea  surface,  and  hydrophone  outputs 
were  recorded  on  magnetic  tape.  While  the  prime  purpose  of  this  experiment 
was  to  measure  ambient  noise  only,  several  "cuts"  were  made  with  a target 
submarine  in  the  test  range  vicinity.  A total  of  eleven  5-minute  test 
events  were  recorded  when  no  traffic  was  present;  two  events  were  recorded 
with  a submarine  on  the  surface  6000  yd  from  the  test  sight,  and  five 
events  were  recorded  when  the  submarine  was  making  a submerged  run  through 
the  test  area.  For  the  present  effort,  however,  analysis  was  performed 
only  on  the  ambient  sea  noise  data. 


' . i.  .<  • 


i&if' 


Data  was  acquired  by  analog  recording  all  unprocessed  hydrophone  outputs 
on  a 14-channel  tape  recorder  with  time  code  and  voice  notation.  Refer- 
ring to  the  equipment  depicted  in  figure  8,  a tape  recording  speed  of  30 
ips  yielded  about  15  minutes  of  test  data  per  reel,  with  the  20  dB  pads 
providing  the  signal  reduction  necessary  for  tape  recorder  compatibility 
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Figure  7 


Gradient  Hydrophone  Calibration  Beam  Pattern 


Figure  8 . Data  Acquisition  Equipment  Setup 


To  ensure  that  a proper  reference  to  input  sound  pressure  level  was  ob- 
tained, the  system  was  calibrated  and  tape  recorded  at  1/3  octave  center 
frequencies.  A minimum  frequency  range  of  150  Hz  to  3 kHz  was  recorded. 

Since  there  was  a possibility  of  recording  radiated  noise  from  targets  of 
opportunity  along  with  the  background  noise,  the  data  were  verified  and 
edited.  As  shown  in  figure  9 , data  recordings  of  an  omnihydrophone  were 


Figure  9.  Data  Editing 
11 


processed  through  an  FQQ  analyzer.  The  resultant  frequency  versus  time 
presentation,  coupled  with  the  time  code  information,  was  used  to  eval- 
uate which  events  truly  depict  background  noise  or  a target  of  opportunity 
To  verify  that  the  spectral  levels  were  correct,  narrowband  frequency 
spectra  were  produced  for  selected  events. 

The  noise  data  were  reduced  into  usable  form  by  first  passing  them  through 
an  analog  filter.  The  analog  bandpass  filters  were  General  Radio  fre- 
quency-selective 8%  filters  (3  dB  bandwidth  is  8%  of  selected  center  fre- 
quency) . For  convenience,  a tap  was  made  of  the  bandpass-filtered  noise 
waveform.  These  noise  data  were  then  converted  into  a suitable  digital 
format  compatible  with  the  storage  capabilities  of  the  computer.  These 
stowed  noise  data  were  then  available  for  use  in  the  simulation  stage  of 
the  investigation. 

SYSTEM  SIMULATION 

The  raw  hydrophone  noise  information  was  converted  into  CTASS  array  per- 
formance evaluation  data  by  simulation.  In  general,  a synthesized  signal 
and  interference  (after  suitable  weighting  dependent  upon  the  detection 
operation  being  simulated  , and  the  hydrophone's  directivity  response) 
were  operated  upon  and  considered  to  be  that  hydrophone's  output  response 
as  if  it  were  receiving  actual  target  signal (s)  and/or  interferences,  plus 
the  background  noise.  This  combining  of  in  situ  noise,  plus  the  synthe- 
sized signal  (s)  would  allow  for  the  evaluation  of  a wide  range  of  signal, 
noise,  array  parameters,  and  geometry.  Data  processing  after  that  point 
paralleled  that  which  would  occur  in  an  operational  system  for  an  individ- 
ual frequency  bin.  Figure  10  shows  the  basic  analysis  process. 

The  first  step  in  the  processing  procedure  was  to  pass  the  hydrophone  out- 
put through  a narrow  band  digital  filter.  The  filtered  output  was  then 
multiplicatively  combined  with  that  from  companion  hydrophones,  and 
averaged  to  develop  the  directional  character  of  the  array.  The  output  of 
the  averager  could  then  be  analyzed  to  form  the  desired  statistical  quan- 
tities needed  to  predict  array  performance.  A detailed  description  of  the 
signal  processing  used  in  the  simulation  program  appears  in  appendix  A. 


In  each  simulation  run,  the  quantity  of  interest  is  the  output  of  the 
integrator  r(T)  (which  depends  on  integration  time  T)  and  its  statis- 
tics (mean  value  and  variance) . Plane-wave  propagation  is  assumed  for 
signal  and  interference,  with  direction  angles  eg  and  0 , respectively. 
Pure  sine  waves  are  assumed  with  peak  amplitudes  Pc  and  frequencies  f 
and  f (1  + AfQ) , respectively.  Initial  phase  angles  are  taken  to  be 
zero  for  convenience.  We  thus  have: 

2 = ^oo^  ' (yl  y-l>  ’ (y2  y-2}  ( > (1) 


m = 4 (1  section)  etc 


n = 2 (BPU) 


= Z (kAt) 
m 

and: 

- k\ml 

where:  T = nAt 

At  = A/D  conversion  sampling  time 
y^  (kAt)  = Noise  + Signal  at  output  of  digital  filter 

The  statistics  are  formed  by  averaging  over  several  T-second  intervals. 
Letting  r^  be  the  output  of  the  integrator  for  a given  array  when  n^ 
samples  are  integrated,  we  define  the  mean  value  as: 


i=l 

and  the  mean-squared  value  as: 


(4) 
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The  variance  is  given  by: 


= r2  - Q2 


We  distinguish  between  the  following  three  cases: 

Noise  only  present  (N) . 

Signal  and  noise  present  (S  + N) . 

Signal  plus  interference  plus  noise  present  (S  + I + N) 
The  output  signal-to-noise  ratio  is  defined  as: 
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(n)0 
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= 1^1 
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where  Ar  is  the  shift  in  the  mean  value  caused  by  the  presence  of  the 
signal : 

Ar  = r„  - r 
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2.  Output  S/N  ratio,  defined  as 
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S+N 


where  o|+N 


is  the  variance  produced  by  signal  plus  noise,  as  given 
by  equation  (5) . 


Input  signal-to-noise  ratio  is  defined  as  the  ratio  of  signal  power, 

2 

(Pg/2 ) , to  the  noise  power  that  is  measured  by  the  omnidirectional  hy- 
drophone. An  appropriate  scaling  factor  must  be  included  to  account  for 
the  effects  of  the  hydrophone  sensitivity,  tape  recorder,  analog  filter, 
A/D  conversion  and  digital  filter  bandwidth.  Input  S/N  ratio  is  referred 

to  as  (li  • 

It  was  anticipated  that  the  following  four  series  of  tests  would  be 
conducted: 

1.  The  bean  pattern,  defined  as  r(©s),  which  depends  on  frequency  f 
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3.  The  interference  effect,  evaluated  first  with  respect  to  the  change  in 
mean  value  caused  by  the  introduction  of  an  interference;  i.e.,  we  will 
calculate : 


Arl  rS+I+N 


- r 


S+N 


(R) 


and,  then,  with  respect  to  the  change  in  (S/N) 


S')  (rS+I+N 

— I = j 


- y (; 


rS+N 

T 


'W 


(9) 


S+I+N 


S+N 


4.  Resolution  was  to  be  made  for  the  special  case  of  equal  signal  and 

interference  amplitude,  Pg  = P^.  For  a given  run,  the  angular  separation 
between  interference  and  target  remains  fixed  (A0  = 0 - 0 = constant) , 

b 1 

but  the  output  is  evaluated  for  several  values  of  0 . This  is  equivalent 
to  rotating  the  array,  or  to  a preformed  beam  system  in  which  multiple 
beam  outputs  are  compared.  We  will  evaluate  Ar  and  as  a function 

of  0 , where  ° 


Ar  = r 


S+I+N 


rN  ' 4rs+I 


and 


ft)  ■ 


( ArS+I)  2 
°2S+I+N 


(10) 


(11) 


Table  1 is  a compilation  of  the  input  and  output  parameters  that  were 
to  be  varied  in  the  tests.  For  each  test,  the  variation  of  the  output 
parameter,  as  affected  by  the  number  of  hydrophone  sections  processed, 
was  to  be  ascertained  for  a fixed  set  of  input  parameters. 
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Table  1 


Test  Summary  (Anticipated) 


Type  of  Test 
Beam  Patterns 

Output  (S/N)q 


Interference 

effects 


Resolution 


Output  Parameters 
r (0) 


i 

Input  Parameters 

2 frequencies 
2 inputs  (S/N) 

0g  from  0°  to  180° 

] 

2 frequencies 
Several  inputs  (S/N) 


2 integration  times 

1 frequency 
1 input  (S/N) 


Several 

(I/S) 

Several 

9I 

Several 

Afo 

1 frequency 

1 input 

2 Afo 

(S/N) 

Several 

A0 

Several 

es 

S IMULATION  RESULTS 


A number  of  CTASS  simulation  runs  were  made  using  recorded  sea  noise  that 
had  been  passed  through  8%  General  Radio  analog  filters  and  converted  from 
analog  to  digital  format,  as  described  in  appendix  A.  Frequencies  of  525 
Hz  and  615  Hz  were  selected.  These  frequencies  were  chosen  after  exami- 
ning the  spectral  plots  of  recorded  data.  It  is  believed  that  essentially 
only  isotropic  sea  noise  was  present  in  the  excerpts  taken  from  the  tape 
recorder. 

Although  the  simulation  program  was  developed  to  handle  multiple  targets 
in  order  to  evaluate  the  resolution  and  interference  discrimination  capa- 
bilities of  CTASS,  the  first  order  of  business  was  planned  to  be  a realis- 
tic evaluation  of  the  output  signal-to-noise  ratio  (S/N)Q.  If  CTASS  did 
not  show  promise  on  the  basis  of  that  criterion,  no  further  simulation 
work  was  to  be  performed.  The  object  of  the  first  set  of  runs  was  to 
evaluate  the  (S/N)  for  the  CTASS  array  consisting  of  various  sections, 
for  "on-target"  conditions;  i.e.,  for  a simulated  target  lying  on  the 
acoustic  axis  or  main-response  axis  (MRA) . Several  values  of  input  sig- 
nal-to-noise  ratios  (S/N)^  were  used. 

The  results  were  very  discouraging.  Contrary  to  what  had  been  hoped,  it 
was  found  that  (S/N)q  decreases  as  the  size  of  the  CTASS  array  increases. 

In  other  words,  for  a given  set  of  fixed  inputs  to  the  simulation  program 
(such  as  sampling  frequency,  signal  amplitude,  etc) , the  output  S/N  ratio 
is  largest  for  M = 2 (basic  processing  unit) , then  decreases  successively 
for  M = 4,  6,  8,  and  10.  This  was  found  to  be  the  case  fairly  consistently 
for  several  values  of  (S/N)^. 

Earlier  analytical  work  done  with  broadband  signal  and  noise  inputs  sug- 
gested that  this  decrease  in  (S/N)Q  might  be  encountered  with  increasing 
number  of  array  elements  ^ ^ . The  mathematical  complexities  did  not  per- 
mit a complete  analysis  of  the  general  CTASS  array.  Only  the  basic  pro- 
cessing unit  (BPU,  M = 2)  was  analyzed  fully.  A comparison  of  the  output 
S/N  ratio  between  M = 2 and  M = 4 was  made  in  reference  (1),  but  a dif- 
ferent definition  was  used  for  output  S/N  ratio.  The  variance  in  the 
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denominator  of  the  S/N  ratio  equation  was  taken  as  being  caused  by 
noise  only,  and  the  fluctuations  caused  by  the  signal  were  disregarded. 
This  was  done  for  two  reasons:  to  make  the  analysis  more  tractable,  and 

because  the  case  of  interest  in  broadband  applications  is  that  of  small 
input  S/N  ratios.  For  very  low  values  of  (S/N)^,  the  variance  with  signal 
present  is  practically  the  same  as  the  variance  from  noise  alone.  The 
output  S/N  ratio  thus  defined  was  found  to  be  equal  to  28.8  (S/N)^  BfT 
for  M = 2,  and  (23.7 (S/N)^  + 96  (S/N)3  + 97  (S/N)4)  BfT  for  M = 4.  where 
Bf  is  the  filter  bandwidth  and  T the  integration  time  of  the  averager. 

The  fact  that  the  (S/N)2.  coefficient  for  M = 4 is  slightly  smaller  than 
the  corresponding  coefficient  for  M = 2 gave  rise  to  the  suspicion  that 
output  SNR  might  decrease  as  M increases,  but  it  was  thought  that  the 
other  terms  present  in  the  M = 4 expression  (which  would  also  be  present 
for  higher  order  systems)  might  cause  the  output  S/N  ratio  to  increase  as 
M increases,  particularly  for  larger  input  S/N  values  (as  are  present  in 
narrowband  applications) . 

For  this  reason,  it  was  proposed  for  the  simulation  program  to  evaluate 
the  output  S/N  first  (which  now  has  to  include  the  signal  fluctuations 
because  of  the  larger  input  signal-to-noise  ratios  used  with  narrowband 
signals)  to  see  whether  it  increases  with  M,  before  any  other  aspects  of 
CTASS  performance  were  to  be  evaluated.  Further  work  was  to  cease  if 
the  early  results  turned  out  to  be  disappointing. 

Figure  11  shows  the  ratio  of  output  S/N  ratios  relative  to  the  basic  pro- 
cessing unit;  i.e.,  we  plot  R vs  M,  where: 

(S/N)  (M  elements) 

(S/N)  o (M  = 2) 

The  input  S/N  ratios  and  the  number  of  samples  (k)  used  in  the  smoothing 
algorithm  are  indicated  by  the  arrows  in  the  figure  for  each  curve.  It 
can  be  seen  that  R > 1 for  M = 4 in  two  cases.  This  is  probably  because 
of  the  finite  averaging  time.  Except  for  these  aberrations,  the  curves 
are  seen  to  exhibit  the  same  trend,  decreasing  with  increasing  values  of 
M.  The  slope  is  particularly  steep  for  small  (S/N) ^ values.  Howover, 
even  for  the  extremely  large  values  ( (S/N)^  = 2 x 10®  and  K = 104),  R is 
seen  to  decrease  with  M. 
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(S/N)  (M  ELEMENTS) 


; 


In  the  noise-free  case,  this  becomes  merely: 


(-)  = 


(rj 


r2"-  (?) 


(5) 


where  ? and  r2  are  the  average  values  of  equations  (2)  and  (3) , re- 
spectively. Since  the  periodic  components  do  not  contribute  anything 
to  the  mean  values,  we  obtain: 


r = P2/2 


(6) 


_ p2  2 / 

^ - <r->  l1  + 


a2  + b2 


But  a|  + k>2  = 2 (l-cos2wT) 


4w2T 


8oi2T2  ' 

Substituting  into  equation  (5),  we  get: 


(7) 


,S.  - ^ '-I-'  /ON 

N 1 - cos2wT 
o 

Note  that  (S/N)  o -*■  00  as  cos2toT  -*•  1.  This  happens  when  T is  very 
close  to  an  integral  number  of  periods  of  the  signal. 

Because  of  round-off  errors,  it  is  not  possible  in  a digital  system  to 
make  T (which  equals  KAt)  exactly  equal  to  an  integral  number  of 
periods  of  the  frequency  before  multiplication.  Hence,  there  will  be 
some  fluctuations  in  the  output  value  caused  by  the  signal,  and  (S/N)q 
will  be  finite  even  in  the  noise-free  case.  In  practical  applications 
this  is  even  more  true,  since  the  frequency  is  never  known  exactly, 
and  since  any  tone  will  have  a certain  amount  of  modulation. 

If  we  write  T as: 

(9) 


T = nT  + A 
o 


(n  = 1 , 2 , ) 


1 2ir 

Where  T = -e—  = — , we  have: 
of  0) 

o 

2u>T  = 2w  ^n  —■  + = 4nir  + 2wA 

cos2wT  = cos(4nir  + 2oiA)  = cos  2wA 

For  2wA  <<  1,  we  can  use  the  approximation: 

_ _ „ „ , (2u)A)  2 

cos  2wA  = 1 - - — = 


(10) 

(ID 

(12) 


22 


So  that : 


(S)  4^ 


2T‘ 


2to2  A‘ 


( 


This  covers  the  case  M = 2.  Equation  (8)  represents  the  output  S/N 
ratio  for  any  value  of  T,  and  equation  (13)  is  valid  for  2u>A  <<  1. 


To  analyze  the  case  M = 4,  we  proceed  in  like  fashion.  We  have: 

t 

T ~r  I cos40JT  dx  c 


r(t)  = i A2 P 4 j cos4tox  dx 


t-T 


where  A^  is  the  angular  sensitivity  of  the  two  additional  gradient 
hydrophones  (tilted  at  +10°  in  our  case,  making  A-^  = cos  10°)  to  a 
target  on  the  MRA.  We  now  obtain: 


r (t ) = 


A1  p4 


toT 


3 , sin  2x  , sin  4x 

8 4 32 


tot 


C 


to  (t-T) 

Carrying  out  the  calculations  as  in  the  previous  case,  we  get: 

3 


(?)2  = (f  A2P4) 


C 


ll  A2p4)2 

1 ♦ < + ^ 

. + bMl 

\ 8 A1  ) 

9(o2T2  ^ 2 

128  ' 

where  a_  and  b_  are  the  same  as  before,  and  a.  = 1 - cos4coT,  b,  = 
2 2 4 4 


sin4toT.  Simplifying  (17)  and  subtracting  (16),  we  obtain: 

2 


9to2T2  ' 


cos2toT  + 


1 - cos4toT 


64 


&) 


(] 


Finally, 


(£,  = 
No 


ito2T2 


1 - cos2<oT  + 


1 - cos4iot 


(1 


Again,  we  see  that  (S/N)q  -*•  00  as  u)T  -*■  mr  (n  — 1,  2,  ...)• 

Letting  T = nTQ  + A as  before  (equation  9) , we  obtain: 

,S.  18  T^ 

' 17  A2 

Thus  we  see  that  (S/N)  has  decreased  when  going  from  M = 2 to 
M = 4.  Carrying  one  step  further,  we  analyze  the  6-element  array 
(M  = 6).  Now  we  get  for  the  output  of  the  integrator: 


(20) 


A?A2P6  wt 


r (t)  = 


12 


/ 


cos6x  dx 


oo(t-T) 

Carrying  out  the  mathematical  details,  we  obtain: 


(21) 


r = i-5-  A2A2P6 
r 32  a1a2f 


(22) 


CM 

| 

fl  + 1 

oj2T2  * 

(1  - cos4wT)  + (1  - cos6wT 

The  output  S/N  ratio  becomes: 


(£,  = 
N o 


400  . . 2 m 2 

225  u T 


_ _ . 1 - cos4wT  . 1 - cosGujT 

1 - cos2u)T  + + : 

25  (45) 2 


(23) 


(24) 


In  terms  of  A,  we  obtain: 


,S.  _ 100  T 
¥o=  131  A 


(25) 
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^ •>  7.  iT* 


We  see  that  the  output  S/N  has  sufferer,  another  decrease  when  going 
from  M = 4 to  M = 6.  If  we  divide  the  output  S/N  ratio  obtained  for 
M = 4 and  M = 6 by  that  for  M = 2,  we  obtain  for  the  ratio  R = 0.53 
for  M = 4 and  R = 0.38  for  M = 6.  These  values  are  indicated  in  fig- 
ure 11  as  asterisks  to  show  the  comparison  with  experimental  values. 

There  is  no  reason  to  suppose  that  this  trend  of  decreasing  output  S/N 
ratio  with  increasing  array  size  would  be  reversed  for  M > 6. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  agreement  between  the  theoretical  and  experimental  results  shown 
here,  coupled  with  the  analysis  performed  on  broadband  signals  (ref- 
erence 1)  which  implied  a decreasing  output  S/N  ratio  as  M increases, 
leads  us  to  conclude  that  CTASS  in  its  presently  conceived  form  is  not 
a promising  array  processing  system. 

It  is  entirely  possible,  however,  that  some  other  way  of  combining  the 
outputs  of  the  gradient  hydrophones  (e.g.,  pair-wise  multiplication 
followed  by  addition  and/or  split-beam  processing) , might  lead  to  im- 
proved output  performance  while  retaining  some  of  the  other  CTASS  ad- 
vantages (high  maneuverability,  directionality,  etc).  The  S/N  ratio 
evaluation  of  alternate  CTASS  configurations  should  be  made  by  ana- 
lytical methods  (as  done  here  for  the  noise  free  case)  to  the  fullest 
possible  extent.  If  one  or  more  such  configurations  appear  promis- 
ing, a simulation  model  comparable  with  the  one  used  here  should  be 
developed.  The  recorded  at-sea  data  collected  during  this  project 
could  be  used.  The  cost  of  such  an  investigation  would  be  modest  in- 
deed, since  a substantial  portion  of  the  necessary  groundwork  has  al- 
ready been  laid. 
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Appendix  A 

CTASS  SIMULATION  AND  DATA  REDUCTION  PROGRAM 

In  order  to  save  on  computing  and  storage  requirements, 
it  is  desired  to  select  as  low  a sampling  rate  as  possible 
(that  was  one  reason  for  using  the  analog  bandpass  filters) . 

It  is  well  known  that  the  sampling  process  causes  changes  in 
the  spectrum  of  the  original  signal.  In  particular,  "aliases" 
are  introduced;  that  is,  other  components  appear  in  the 
spectrum  after  sampling  that  were  not  present  before  sampling. 
The  mathematical  relationship  between  the  spectrum  of  the 
sampled  signal  Fs(jw)  and  the  original  signal  F ( j oj ) is 

oo 

F (jw)  = I ak  F ( jw-kojg)  (1) 

k=-®> 

where  w = 2vf  (f  = sampling  frequency  = 1/At) , and  a, 
depends  on  the  type  of  sampling  (for  impulse  sampling 
assumed  here,  a^  = constant) . 

Figure  1 illustrates  the  sampled  spectrum  for  a narrow- 

band  signal  centered  at  f and  a sampling  frequency  of 

f = % f . The  shaded  portions  correspond  to  the  spectrum 
s 9 o 

of  the  analog  signal  (before  sampling) . 


Nyquist 

Interval 


Figure  1.  Spectrum  of  Sampled  Signal 


A-l 


It  can  be  seen  that  the  sampled  signal  has  a rather  compli- 
cated spectrum.  Fortunately,  i.ne  only  portion  of  the  spectrum 
which  is  of  direct  importance  is  the  part  lying  inside  the 
Nyquist  interval 


For  our  example,  this  is  the  "second  alias,"  namely,  the 
positive  frequency  portion  shifted  by  -2fg  and  the  negative 
portion  shifted  by  + 2fg. 


It  can  be  seen  that  we  have  thus  effectively  "band- 

shifted"  the  data  down  to  a lower  frequency  by  sampling  the 

data  at  a rate  lower  than  the  frequencies  contained  in  the 

original  signal.  In  order  to  minimize  distortion  between 

aliases,  it  is  desired  to  have  the  distance  between  them  as 

large  as  possible.  At  the  same  time,  we  wish  to  keep  our 

4 

sampling  rate  low.  The  particular  choice  f = ^ fQ  is  a 
good  compromise.  The  interferences  caused  by  the  nearest 
alias  on  either  side  of  the  spectral  component  within  the 
Nyquist  interval  are  made  equal  (for  a symmetrical  spectrum, 
which  is  assumed  here)  by  this  choice  of  sampling  frequency. 

Since  the  bandwidth  B of  the  analog  filters  (hence  of 

the  noise  spectrum)  is  .08  f , and  since  the  distance  between 

o 

aliases  is  2fo/9  = 0.2222  •••  fQ,  the  nearest  alias  is 
approximately  5 half-bandwidths  away  from  the  spectral  com- 
ponent of  interest.  This  corresponds  to  a distortion  due  to 
aliasing  of  less  than  20  db  in  amplitude,  and  can  be 
considered  as  negligible. 


In  order  to  achieve  locally  high  input  signal-to-noise 
ratios  even  for  low-level  signals,  it  is  desired  to  use 
narrowband  filters.  Because  of  the  difficulty  in  designing 
very  narrow  bandpass  filters,  it  is  more  convenient  to  band- 
shift  the  spectrum  down  to  base-band  and  use  a low-pass 
filter  to  obtain  the  same  end  result.  This  is  done  by  split- 
ting the  signals  up  into  quadrature  components,  as  will  be 
shown . 

A-2 


. ..  • •-•••'•.vjgrw'rVTTV' 


The  target  and  interference  signals  need  not  be  filtered, 
since  a sinusoidal  signal  of  known  frequency  is  assumed  and 
is  simulated  in  each  case.  In  an  actual  system  with  real 
signals,  the  output  of  the  ith  hydrophone  can  be  written  as 

v^(t)  = n^(t)  + A^(0)P  cos  (u)Q  [t-T^]+a) 

+ A±  ( 4> ) Q cos  [(/(t-cKjt  0]  (2) 

where 

n^ (t)  is  the  ambient  noise 

P is  the  peak  amplitude  of  the  target  signal 
Q is  the  peak  amplitude  of  the  interference 
signal 

wo  is  the  frequency  of  the  target  signal 
u)'  is  the  frequency  of  the  interference  signal 
0 is  the  target  angle  relative  to  the  array  axis 
$ is  the  interference  angle  relative  to  the 
array  axis 

A^  ( ) is  the  angular  beampattern  of  the  ith 
hydrophone 

is  the  geometric  delay  of  the  target  signal 
a.  is  the  geometric  delay  of  the  interference 
signal 

a,  B are  arbitrary  phase  angles 


Since  the  noise  is  assumed  to  be  narrowband  Gaussian  noise, 
it  can  be  written  as 

ni(t)  = n.^  (t)  cos  wot  + ni  (t)  sin  u>ot  (3) 

c s 

Without  loss  of  generality,  the  initial  phase  angle  a can  be 
dropped,  and  Eq.  2 can  be  written  (using  Eq.  3)  as 
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v.  (t)  = (n.  (t)  + P + Q ) cos  a)  t 

1 X c c o 

c 

+ (n.  (t)  + P + Q ) sin  w t (4) 

s 

= v.  (t)  cos  a)  t + v.  (t)  sin  w t 
l o l o 


This  is  now  in  the  form  of  the  two  quadrature  components, 

(t)  cos  wQt  and  (t)  sin  u)Qt.  It  can  readily  be  shown 
c s 

that  when  v^(t)  is  multiplied  by  2 cos  u>Qt  and  the  output 

is  passed  through  a low-pass  filter,  the  filtered  version 

of  v,  (t)  is  obtained.  The  filtered  version  of  Vi.  (t)  can 
c s 

similarly  be  gotten,  and  the  filtered  form  of  v^(t)  can  be 

constructed  by  multiplying  vic(t)  by  cos  wHt,  vig(t)  by 

sin  u)..t  is  some  convenient  frequency,  not  necessarily 

n n 

equal  to  u)Q)  and  adding  the  two  signals. 

This  procedure  will  be  used  in  our  analysis  program, 
except  that  the  frequency  f is  used  instead  of  fQ,  and 
that  the  target  and  interference  signals  will  be  simulated 
as  levels  Pc,  Qc,  Pg , and  Qg , rather  than  as  sinusoidal 
functions.  This  simplifies  the  simulation,  but  is  equiva- 
lent. In  fact,  it  also  eliminates  any  phase  errors  which 
might  result  from  inaccurate  knowledge  of  the  center  fre- 
quency fc-  The  diagram  of  Fig.  2 illustrates  the  procedure. 

It  was  drawn  for  an  arbitrary  channel.  Naturally, 
for  the  basic  processing  unit  the  outputs  of  the  omniphone 
and  the  gradient  phase  are  combined  to  replace  the  omni 
channel.  As  many  as  10  outputs  are  multiplied  at  the  far 
right  of  Fig.  2 . Because  this  gives  rise  to  frequencies  up 
to  10  o>H,  care  must  be  taken  that  the  highest  frequency  is 
still  within  the  Nyquist  band 


Block  Diagram  of  Digital  Analysis  Program 


i.e.,  we  require 

10  fH  + r-  s 5s  <5> 

where  B'  is  the  bandwidth  of  the  digital  filter.  Since 
there  will  also  be  frequencies  of  2f„,  4f„,  6f„,  and  8f„ 

rl  rl  n n 

present,  we  also  require  that  there  be  sufficient  separation 
between  them,  or  that 


B’  < f. 


(6) 


The  following  limits  satisfy  these  requirements: 


f H * 50  and  B'  - 100 


This  gives  (using  the  equality) , 


B'  = £„/ 2 < £„ 


and,  since  fs  = 5 fc 


10fH + r-  - r + iM  < I £o 


An  input  frequency  of  f = 500  Hz  thus  allows  us  to  go  up 
to  a bandwidth  of  B'  = 5 Hz. 

The  notation  used  is  listed  below  and  a complete  set 
of  equations  for  the  simulation  program  is  presented.  The 
constants  are  used  to  neutralize  the  effect  of  different 
hydrophone  sensitivities. 

n^Q(kAt)  = digitized  noise  from  omni  phone 
digitized  noise  from  gradient  phone  with  n = 0 
digitized  noise  from  gradient  phone  with  n = +10° 


n (k)  = 
oo 


noo(k)  = 


n1(k)  = 


. , f - 


A- 6 


n_1 (k) 

— 

digitized 

noise 

from 

gradient 

phone 

with 

n 

= 

-10 

n2  (k) 

= 

digitized 

noise 

from 

gradient 

phone 

with 

n 

= 

+ 20 

n_2 (k) 

= 

digitized 

noise 

from 

gradient 

phone 

with 

n 

= 

-20 

n3  (k) 

= 

digitized 

noise 

rom 

gradient 

phone 

with 

n 

= 

+ 30 

n_3 (k) 

= 

digitized 

noise 

from 

gradient 

phone. with 

n 

= 

-30 

n4  (k) 

= 

digitized 

noise 

from 

gradient 

phone 

with 

n 

= 

+ 40 

n_4 (k) 

= 

digitized 

noise 

from 

gradient 

phone 

with 

n 

= 

-40 

k 

= 

kAt;  At  = 

1 

r~  ~ < 

s 

9 

if 

o 

2.25 

o 

= speed  of  sound 


d. 

— - sin  0 
c 


u • 

= — sin  <p 


i = ±1  to  ±4 


-nB'At 


a2  = 1 - 3 


n (k) 


[K;o^o{k)  + Koonoo(k)] 


no  (k) 


nQ(k)  cos  tock 


dQ  (k)  = (k-1)  + a2no  ^ 

o o c 

e (k)  = a, e (k-1)  + a2d  (k) 


Digital  Filter  (I) 


xo  (k) 


£ (k)  + 5(l+cos0)  + §(l+cos  4>)  cos  (Aw  k + 8) 

°c  2 2 ° 


y0  no 


x (k)  cos  to„k 

O rl 

C 


...  '.V-..-?  ; •*. 


noo  (k)  = 2n(k)ooKoo  coswck 
c 

d , Jl  . , _ . 

oo  oo  as  in  (I) 
c c 

x (k)  = H (k)  + P COS0  + Q cos<f>  cos(Acj  k + 8) 
oo  oo  o 

c c 

y (k)  = (k)  cos  u>„k 

J OO  OO  H 

c c 

yQo(k)  = yoo  (k)  + yooJk) 
c s 

noo  (k)  = 2noo(k)Koo  sin  wck 
s 

d , l as  in  (I) 
oo„  oo 
s s 

x (k)  = l (k)  + Q cosij)  sin(A(i)  k + 8) 
oo  oo_  o 

s s 

yoo  (k)  = xoo  (k)  sin  wHk 
s s 


A- 8 


ni(k) 


(i  = ±1,  ±2,  ±3,  ±4) 


n.  (k) 


[K.n.  (k)  cos  co  k]  x 2 
IX  c 


n . 

1 s 

(k) 

= [Kini(k)  sin 

a>ck]  x 2 

d. 

(k) 

- ii  (k) 

as  in 

(I) 

xc 

c 

d. 

Xs 

(k) 

-*•  (k) 

s 

as  in 

(I) 

X . 
1 

(k) 

n 

H- 

l 

XT 

+ P cos  (0-n^)coscooT^ 

+ Q cos$  -n^)  cos  (Aojok  + wQa^) 


x.  (k)  = (k)  + P cos  (0-ni)  sina>oTi 

1s  s 


+ Q cos  ) sin  (AcoQk  + coot^) 


y i (k)  = xi  (k)  cos  coHk 
c c 

y.  (k)  = x.  (k)  sincu  k 
xs  s 

yi (k)  = yi  (k)  + yi  (k) 
c s 


Now  form 


zo (k)  = yo(k)yoo(k) 

(BPU) 

M = 2 

Zi()c)  = z0(k)y1(k)y_1(k) 

(M=4 ) 

- 1 section 

z2(k)  = z1 (k)y2 (k)y_2 (k) 

(M=6) 

- 2 section 

z3(k)  = z2(k)y3(k)y_3(k) 

(M=8 ) 

- 3 section 

z4(k)  = z3 (k)y4 (k)y_4 (k) 

(M-10) 

- 4 section 

1 L 

r.  (L)  = f l Z-iOO 
3 L k=l  3 

print  running  average 
and  store  for  L = L 

n 

«■* . ' . . 77*  ■ 

- ' . J • - -V 

<+•  •> 

* ‘ .V  **■  1 


t 


If  we  integrate  over  two  periods  of  the  lowest  frequency 
of  interest  (2  co^)  , we  let 


2f  H f H 


If  f H 50  and 

f = 
o 

‘i  r:  ..  1 

9 f o r At  ■ f 

y 

4f 

f 

s 

o 

o 

rn  _ 1 — 

50 

L - X - 

50/fo  _ 

200 

■ f H 

f ' 
o 

n At 

! fo 

9 

Should  get  very  little  fluctuation  after  an  average  of 
about  10  to  20  samples. 

Noise  Reference 


(k)  = 2n  (k)  K ’ 

cos  co  k 

OO  oo 

c 

(k)  - aldoo  (k-1} 
c 

+ a2noo (k)  1 

(k)  = al£oo  (k_1) 

+ a2doo(k)  J 

c 

c 

- £oo  <k) 

Dig.  Filter  (I) 


y ™ (k)  = coswHk 


noos(k)  = 2noo(k)K;o  sin  uck 

doo  (k)  = aldoo  (k_1)'  etc- 
s s 

loo  (k) 
c 

. I 

X = 0 

OO  OO  r 

s c 

yoo(k)  = Xoo  (k)  sinuHk 
s s 


Dig.  Filter  (I) 


y00(k)  = y«o  + y™  (k) 


- tyoo<k>> 


A-10 


■ **»r- 


The  first  set  of  equations  with  the  subscripts  o and  oo 
apply  to  the  basic  processing  unit.  The  next  set  of 
equations  are  for  added  sections.  Because  many  of  the 
equations  are  repeated  (only  the  subscripts  differ) , only 
the  general  subscript  i was  used.  The  digital  filter 
equation  denoted  by  (I)  is  repeated  again  and  again  and 
is  thus  subsequently  merely  referenced  by  (I) . 
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